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ABSTRACT 

Thermodynamics, when formulated as a dynamic action principle, allows a simple 
and effective integration into the General Theory of Gravitation. 



I. Introduction 

Extensions of Einstein's theory of gravitation were sought, and found, from 1916 
onwards. This paper deals with the incorporation of thermodynamics into General 
Relativity, the relativistic extension of thermodynamics. We begin by declaring the 
aims and the limitations of this work. 

The prototype of a unified field theory is the unification of the theories of Maxwell, 
Einstein and Dirac that is defined by the lagrangian. But 'electrodynamics' is a vast 
edifice that includes, besides the elementary interactions of photons and electrons, the 
whole field of electromagnetic interactions of matter, including most of solid state 
physics. The Maxwell-Einstein system does not encompass electromagnetism in this 
widest sense; it is principally a relativistic theory of elementary particles, and it is easy 
to understand why: Both theories are formulated in terms of action principles, all it 
takes is to add the actions of both theories. It is significant that all unified theories 
studied so far are based on action principles, and this even applies to some propos- 
als for combining Thermodynamics with General Relativity (Taub 1954, Bardeen 1970, 
Schultz 1970). Modern astrophysics is not based on an action principle, and yet Ein- 
stein's action has been invoked many times, since Eddington's book of 1926. The theory 
presented here has been used (Fronsdal and Wilcox 2011a, 2011b) to model the dark 
matter distribution of the Galaxy. 

It is evident that a unification of thermodynamics with General Relativity must 
begin by formulating thermodynamics in terms of an action principle. To lay the foun- 
dations of thermodynamics in terms of an action principle is anything but new, for it 
was the inspiration and guiding light of Gibbs' famous 300 page paper (Gibbs 1875). 

The aim of this work is to give a formulation of thermodynamics in the framework of 
eulerian field theory, where the dynamical variables are densities, pressure, temperature 
and entropy, all fields over space and time, that allows for its integration into General 
Relativity. The solution proposed in this paper is restricted in some ways; applications 
include only adiabatic processes. It is a start. 
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The potential benefits of formulating a fully fledged action principle for thermo- 
dynamics extend far beyond the present context; perhaps the principal advantage is 
that each successful application contributes constraints on the lagrangian that must be 
respected in subsequent applications. One should aim at a "lagrangian for hydrogen" 
and a "lagrangian for water" that would encompass all properties of those substances, 
in all their phases. 

In the strictest sense thermodynamics is the study of equilibria, but of course we 
shall not be confined to that narrow context. In the first place, to make meaningful 
contact with General Relativity, we have to deal with localized thermodynamics, some 
times called Extended Thermodynamics, and Hydrodynamics. In the second place we 
shall have to deal with systems with several components and several phases, and with 
a general class of equations of state. 

It is clear that the focus on an action principle limits the scope of the theory. 
The most severe limitation of theories based on a variational principle is that the flows 
have to be at least piecewise irrotational and determined locally by a velocity potential, 
although limited vorticity can be accommodated in the manner of Onsager's theory 
(Onsager 1949).* A specific lagrangian determines an adiabat and applies exclusively 
to adiabatic fluctuations. That this is useful is due to the fact that, as a rule, non- 
adiabatic changes take place on a very long time scale. Most treatments of continuous, 
non adiabatic changes assume that the time development proceeds through a sequence 
of adiabatic equilibria. In the case of the simplest systems the adiabats may be indexed 
by a single parameter, related to the entropy. A single expression, containing a free 
parameter, then applies to a family of systems interrelated by the exchange of heat. 
Recognition of this fact, together with the approximate additivity of lagrangians of 
composite systems, is a promising starting point for a study of the entropy and the non 
adiabatic processes of heterogeneous systems. That is; the action principle provides 
a convenient framework within which to study entropy, the most difficult concept in 
thermodynamics. (See Fronsdal 2011, Fronsdal and Parthak, 2011.) 

1.1. Dynamical variables and thermodynamical potentials 

The dynamical variables include the metric, the temperature, the pressure, the 
entropy, one or more densities, and one or more irrotational vector fields that represent 
flow velocities. 

The laws of classical, equilibrium thermodynamics are expressed in terms of poten- 
tials U, F, H and G (Massieu 1869). Each is defined as a function of specific variables 
(the natural variables), namely 

U(S,V), F(T,V) H(S,P), G(P,T). 

According to Gibbs, a thermodynamic system is completely defined by any one of these 
functions; an explicit expression for a potential in terms of its natural variables is called 
a fundamental relation. The four potentials are related by Legendre transformations 
and the fundamental relations for a given system are equivalent. 

* Points of vortices are omitted from the spacial manifold and the velocity potential 
is a section of a fiber bundle over the manifold. 
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I. 2. The energy for the simplest system 

The most fundamental concept introduced by Gibbs is the representation of the 
states of a system in terms of surfaces in a real space with coordinates V, T, P and S. 
Gibbs based the entire theory on a principle of minimal "energy" and on a subsidiary 
principle of maximum entropy. The energy is subject to independent variations of 2 of 
the variables, the other 2 held fixed at arbitrary values, the four variables are a priori 
independent. The Gibbsean surfaces are 2-dimensional surfaces in four dimensions, 
determined by the variational equations. Physical configurations are points on this 
surface. 

The action is thus a function of the four independent variables. To define the system 
one specifies the potentials and obtains 2 independent relations among the variables that 
in turn define the Gibbsean two dimensional surface. 

The relations obtained by variation of the generic action must have universal va- 
lidity and each of the following definitions of "the energy" satisfies that criterion, 

E 1 = F(V, T) + ST + VP, E 2 = H(P, S) - ST - VP, 

E 3 = G(T, P)+ST- VP, £ 4 = U(V, S)-ST + VP. 

According to Gibbs, in the description of with adiabatic processes the energy is to be 
varied with entropy and pressure held fixed; this suggests identifying the action with 
Ei, then independent variation of V and T give the standard relations 

dF n dF n 

_ + s = , ^ + p = o. 

The system is thus defined by the free energy and the values chosen for P and S. To 
encode the properties of the system in the expression for the free energy F, rather than 
U or G, is often the preferred choice, see for example Rowlinson (1970). It is a natural 
choice that, by fixing the entropy, focuses on adiabatic processes. 



1.3. Localization 

The local extrapolation of thermodynamics seeks to promote these relations to field 
equations that describe spatial but, in the first instance, time independent configura- 
tions. Until further notice we consider systems with one component and a single phase. 
The function F and the variable S are interpreted as specific densities, the variable V 
as a specific volume. The mass density is p = 1/V and densities /, s are defined by * 

f(p,T) = pF(p,T), s = pS. 

* The representation of the entropy density s as a product pS is meant to imply 
that the specific entropy density S is independent of the dynamical variables. This is 
a specialization, to a wide class of systems, that is adopted throughout this paper. For 
a discussion of some simple systems that do not belong to this class see Fronsdal and 
Parthak 2011. 
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The energy density is 



%:= [ h + VP, h = pv 2 /2 + f(p,T) + sT, s = pS, (1) 

with p, S, T and P treated as independent variables. The vector field v is interpreted as 
a flow velocity. Variation of the field T gives 

i+*=0. (2) 

For a one component system, in the case that the specific entropy density is uniform, 
this is the usual polytropic condition. In the case of multi component systems it has a 
significance that by far exceeds what is usually recognized. 

Local variation of p (variations that vanish at the boundary of the domain E), with 
S and P held fixed, gives 

- + P J = (constant). 
p dp 

The local thermodynamic pressure is defined by 

and the last equation reduces to the hydrostatic condition. 

A h + P n 
grad = U. 

P 

A complimentary variation of the density, with the mass fixed but the volume not, gives 
the result that the thermodynamic pressure must coincide with the "external" pressure 
P at the boundary. In the interior the thermodynamic pressure is balanced by external 
forces. 



I. 4. Dynamics 

We are going to extend the action principle to incorporate the hydrodynamic equa- 
tion of motion for the velocity field, as well as the equation of continuity. We begin by 
introducing the Poisson bracket; for two functionals / and g, 



{f,g} = 




Then the two fundamental equations of hydrodynamics take the form 

P = { H , P} = -div(pv) 

and ^ 

v = {H, v} = —grad — . 



4 



The first equation is the equation of continuity and the second equation will be 
related to the Bernoulli equation. The energy functional as defined by Eq.(l) is thus 
a fully fledged hamiltonian, although the temperature does not (yet) have a canonical 
conjugate field. 

We propose going a step further, passing to a lagrangian formulation. As far as we 
know, this can be done only if the velocity field is irrotational, expressible as 

v = — grad <£>. 

The correct form of the kinetic part of the action was discovered by Fetter and Walecka 
(1980). We define an action and a lagrangian density by 

A = J dtd 3 x£ + VP, C = p$-h, h = pv 2 /2 + f(p,T) + sT. (3) 

This makes $ canonically conjugate to the density p. The equations of hydrodynamics 
are now variational equations. Variation of $ and p gives 

dh ■ dh 

p = ~ai>' = ~dp' 

The first is the equation of continuity and the other is the integral form of Bernoulli's 
equation (see below). 

The lagrangian of Fetter and Walecka served only for hydrodynamics, since the 
temperature had been eliminated by means of an equation of state. In Gibbs' thermo- 
dynamics the temperature as an independent dynamical variable. The variation of T 
gives as before 

f +S =0. (4) 

This 'polytropic relation' is to be regarded as a relation satisfied by the dynamical 
variables p and T that acquires substance once the (specific) entropy has been assigned. 
The function / + sT is a function of three independent variables and 

... d(f + sT) Arr d(f + sT) d(f + sT) 
grad (f + sT) = — grad T H — grad p H — grad S. 

The first term vanishes by the equation of motion and the third term vanishes if S is 
uniform, hence 

A d{f + sT) 

pgrad = grad p. 

dp 

Consequently, "Newton's equation" for the acceleration (the Bernoulli equation) is ob- 
tained by variation of the density, is 

Dv 

p— = -grad p. 
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This thermodynamic action principle, based on the hydrodynamic lagrangian of 
Fetter and Walecka, allows for incorporation into the framework of General Relativity, 
as shall be shown in Part II of this paper. That is not the only advantage of this 
formulation of thermodynamics. The following sections give a very short account of 
some other applications and generalizations. 



!.5. 5. Non adiabatic processes 

The lagrangian applies to adiabatic changes only. Relaxation is non adiabatic and 
irreversible. What makes it possible to talk about adiabatic dynamics is the fact (the 
assumption) that relaxation takes place on a much longer time scale. The variational 
equations determine the temperature as a function of the other variables, but there is no 
equation for the time derivative T. Over short time intervals the temperature follows the 
variations of the density, over long time intervals it is determined by the heat equation. 
To go further we should incorporate the heat equation into our framework, which is 
difficult since this equation describes irreversible phenomena and because the process 
involves external agents. 

There is an alternative, and that is to incorporate additional degrees of freedom 
into the theory. For example, cooling due to electromagnetic radiation can be made 
adiabatic by incorporating the electromagnetic fields and the associated dynamics. Un- 
der favorable circumstances the interaction can be incorporated into the entropy term, 
freezing the electromagnetic field results in a fixed value of the entropy. 

Some insight can be obtained on a case by case basis, especially by a study of 
mixtures and phase transitions. 



I. 6. Generalizations 

Here we describea fairly general situation that can be encompassed by the theory 
in the present stage of development. 

Consider a system described by a number of independent densities, a number of 
independent, irrotational vector fields, temperature, entropy and external pressure, and 
confined to a domain E in space. We suppose that we have an expression for the free 
energy density f(pi, P2, •••',T) and propose the following lagrangian density 

n 

C = Y^pjfa - vj 2 /2 - </>) - f( P i, P2, Pn, T)-TJ2 PiSi ~ |t*. (5) 

3 i 

We have included the newtonian potential <p, in recognition of the fact that the p J s 
are interpreted as mass densities, and the Stefan-Boltzmann pressure, reserving com- 
ments for Section III. The first sum is over a set of independent densities, one for each 
independently conserved density. A velocity potential $j is canonically conjugate to 
each conserved densitiy pj. (Fronsdal 2011). The second sum is over all independent 
densities, independently conserved or not. The central problem is the determination of 
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the total free energy density. The most primitive guess and a good approximation in 
some cases is that a first approximation is obtained by summation over the components, 



f(pi, P2, T) = fiiPi, T ) ? 



(6) 



The specific densities Si, Si,... are usually taken to be uniform; that is, spatially 
constants. The appearance of an arbitrary constant in any adiabatic lagrangian is 
expected, for it must be possible to add heat and this must be reflected in the equations 
of motion. But the appearance of n arbitrary constants (one for each component) 
requires that we specify the path that the system will follow in this n-dimensional 
"entropy space" . For any adiabatic system they are determined by the properties of 
the equilibrium configurations; in certain circumstances the Gibbs-Dalton hypothesis 
may serve to fix the ratios, although corrections are required to account for properties 
of real gases. For a system with chemical reactions, as in the case of the molecular 
dissociation of hydrogen, it was shown that the addition of heat to the compound 
system leads along a path along which Si = S2 is fixed. With this assumption the Saha 
equation (Saha 19) may be derived by variation of the the lagrangian with respect to 
the densities. In the case of a saturated van der Waals gas, the path is defined by the 
Maxwell rule of equal areas, derived from the subsidiar principle of maximal entropy, 
with the result that T(S2 — S2) = e, the evaporation energy. In this case too, non 
adiabatic (but reversible) changes are taking place in which the entropy is changing and 
both cases provide some insight into the use and the concept of entropy. 

1.7. 7. Example, hydrogen atmospheres 

Hydrogen is present in most stars, the molecular form at low temperatures and the 
atomic form at higher temperatures. Here we consider the phenomenon of molecular 
dissociation, below the ionization temperature. The hamiltonian density, in a first 
approximation, is 



The index is 1 for the molecular component and 2 for the atomic component, e is 
the molecular binding energy. The entropy is included as Si = —1Zi\nki,i = 1,2. 
Variation of the densities, with p = p\ + pi fixed, leads to the following equation for the 
concentration r = pij ' p at fixed density 



h = p(X + v 2 /2 + </>)+ ep 2 + ^ n iPi ln 



P 



1 — r 



—T- 5 e~ e/nT 
ep 



This is Saha's equation. It was derived under the assumption that Si = S 2 . 
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8. Phase transitions 

For a pure substance some simple phase transitions and the associated dynamics 
is described by using the free energy of a van der Waals gas. As the temperature is 
lowered (non adiabatically) to the domain where the density is no longer determined 
by temperature and pressure the system behaves in some respects like a mixture, but 
surface tension and gravity leads to separation of the phases. The theory describes 
the equilibrium configurations and adiabatic dynamics but not the process whereby the 
system gives up or gains heat from its surroundings, nor the process of separation of the 
phases. In the saturated, separated phase Eq.(6) is valid. Variation of the densities gives 
the pressure and locates the critical point but it does not yield Maxwell's rule. Maxwell's 
rule can be derived from Gibbs complimentary variational principle of maximal entropy. 

The question of the validity of Eq.(6) is vital in other contexts, such as the propa- 
gation of sound and the determination of critical parameters in gaseous mixtures. The 
latter problem is often formulated in terms of the free energy, while sound speeds are 
derived from very different considerations. This brings us to what is without a doubt 
the most radical consequence of adopting the action principle for all applications of ther- 
modynamics: Every successful application gives information about the lagrangian. For 
example, the study of sound propagation in mixtures at normal temperatures suggests 
corrections to the formula (6) by interaction terms of the very specific form (Fronsdal 
2011) 

The best account of critical parameters of mixtures have been obtained by a modification 
of the van der Waals formula for the pressure, but preliminary results based on the 
following formula for the free energy are very encouraging, 

f(pi,p 2 ,...;T) = ^ fi(pi,T) + a tj {piPj) k . (6) 

i 

The particular form of the interacting term is suggested by observation of the excess free 
energy. When more data on sound propagation at low temperatures become available 
it will be imperative to reconcile observations of critical phenomena with sound speeds 
in mixtures, using the same expression for the free energy. 

The ultimate goal of thermodynamic phenomenology should be to determine the 
lagrangian for each substance, including mixtures of all kinds and including all phases, 
and accounting for a large family of adiabatic properties. 
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II. Thermodynamics and General Relativity 

Following Taub (1954), Schultz (1970) and Bardeen (1970), we propose a partial 
unification of thermodynamics and General Relativity. But unlike previous suggestions 
this one is based on, indeed defined by, a specific lagrangian density, 

£tot = g^^, (R + A) + £ m attten 

where the second term is a relativistic extension of the lagrangian density in Eq.(5), * 
\ i n 

^matter = PA^^M^W ~ ^) ~ /(Pl»P2, Pn,T) -T^PiSi + -T 4 . 

3 1 

The radiation term may perhaps be included in the free energy, but we prefer to add it 
explicitly. The first sum extends over a subset of l...n and includes the densities that 
are separately conserved. The constant A is the cosmological constant and in the case of 
an astrophysical system that extends to infinity it is interpreted as the pressure at the 
boundary at infinity. It is a simple extension of the non relativistic lagrangian to which 
it reduces when in the kinetic term one expands 

il>i = c 2 t + 

and takes the limit as c tends to infinity. Of the dynamical metric there remains, to 
this order, only the Newtonian potential 0. Up to terms of order 1/c 2 , , 

<7 MI Vi, M Vv - c 2 = c" 2 (l - 20/c 2 )(^) 2 - (grad^) 2 - c 2 = 2®, - (grad$,) 2 - 20. 

The action is relativistically invariant (invariant under diffeomorphisms) and it 
follows that the Bianchi identities are satisfied by virtue of the equations of motion. 

The matter energy momentum tensor is 

i 

Hence £, on shell, is the pressure, just as it is in the non relativistic theory. The 
expression is similar to that of Tolman (1934), with some differences: 

1. Any number of densities and vector fields can be accommodated. 

2. The charges that are conserved in the non relativistic theory remain conserved 
in General Relativity, 

<9 M J; = 0, J? = y/=g<ril> jtV . 
* An example of this action principle appears in Fronsdal (2007) . 
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3. The lagrangian is determined by the non relativistic free energy and the inter- 
pretation is the same as in the non relativistic theory. Chemical reactions and changes 
of phase can be described precisely as in non relativistic thermodynamics. 

4. All the equations of ordinary thermodynamics and hydrodynamics are recovered 
in the non relativistic limit. 

5. The simplest examples (ideal gas spheres, radiation neglected) lead to equations 
of motion that are almost identical to those used in simple stellar models. 

6. The black body radiation pressure and energy is included by simply adding the 
term.T 4 /3 to the lagrangian density. This is appropriate in a region where radiation 
is believed to be in equilibrium with matter, as in the photosphere of the Sun. It 
raises or lowers the effective adiabatic index towards the asymptotic limit n = 3 as the 
temperature goes to infinity. It is an improvement of the approximation introduced by 
Eddington (1926) who postulated that the ration p gas /Pradiation is uniform. 

7. The temperature is an independent field, fixed by the equation of motion, Eq.(2). 
This is important in the case that radiation is taken into account for this affects the 
polytropic relation. It leads to a correct inclusion of radiation pressure and energy and 
obviates the need to make Eddingtons assumption. 

8. As to the interpretation of the variables it should be noted that the pressure 
is the on-shell value of the lagrangian density, so that the term —pg^ in the energy 
momentum tensor is the same in both theories, but the dynamics is in the expression 
for the lagrangian and it is lost when the lagrangian is replaced by its value on shell. 
As to the other term, it can be compared to Tolman's expression only in the case of a 
simple system. There is no conserved energy density in General Relativity, but in the 
non relativistic limit 

• • • dC 

pi p 1 p = 1 p =h + C, 

dip 

and on shell it is h + p, where h is the hamiltonian. Hence Tolman's p corresponds to 
the hamiltonian density. An essential difference is that Tolman's formula deals with on 
shell quantities that obey no dynamics besides what is implied by the Bianchi identities. 

9. The entropy enters in precisely the same way as in the classical theory. Any 
advance in the conceptual or operational meaning of entropy is transferred to the rela- 
tivistic theory, without additional complications. 

To all this must be added a discussion of the correct identification of mass density. 
The only effective definition of mass in General Relativity is that obtained in any empty 
region of space time, where the Schwarzschild metric can be directly observed. It is 
contained in the form of either goo or g rr . This definition can be usefully extended to 
regions of low matter density, thus for example 

M(x) = r -(c 2 g 00 -l). 

If the metric is determined locally, by the observation of orbital velocities, then this is the 
correct definition, for non relativistic dynamics is affected by g 00 , much less by the radial 
component of the metric. Yet it is customary to define the mass density in terms of g rr , 
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in a manner that is inconsistent with the transformation law of a density. If there is a 
useful definition of mass density then it is the time component of the conserved current. 
There is no conserved current in Tolman's approach; that theory abandons one of the 
key concepts of hydrodynamics. The integrated density is an invariant characteristic of 
a star, but it is not clear that the value of this quantity must agree with the mass as 
determined, locally or at infinity, by the Schwarzschild metric. 



III. Example, dark matter 

An interesting equation of state for stellar interiors, first proposed by Stoner( ) and 
Chandrasekhar( ), is based is based on a degenerate Fermi gas. A variant was used in 
a celebrated paper by Oppenheimer and Volkov( ) and an account of it may be found 
in Landau and Lefschetz( ). Another variant of this equation of state was used recently 
to model the dark matter in the Milky Way. (Fronsdal and Wilcox 1101a, 1101b) 

Let / be the internal energy density. Little or no attention will be given to the 
temperature; it may be assumed, either, that it is constant, or that the entropy is zero; 
in either case / is principally a function of the density. Variation of the action with 
respect to the density gives, in the case of a spherically symmetric metric, 

1 (cV „ -!)=:,= |. 

Taking the gradient leads to the usual hydrostatic conditions relating pgrad^> to the 
gradient of the pressure 

df . 
P = P d-p- f - 

That is, the gravitational potential is identified with the chemical potential. We note 
that 

dp 

Examples studied by Fronsdal and Wilcox (1011a) have (omitting constants) 
p(4>) = sinh 4 (f). The non relativistic approximation leads to the the Laplace equation, 
A0 = p, 

r 2 d r r 2 d r <j) = sinh <j), 
a "sinh-Emden equation" that is solved exactly by 

0(r) = ln(l + 6/r). 

It can be used to account for the observed rotation curve of visible stars in the Milky 
Way. It predicts nearly constant orbital velocities up to the distance r = b. The solution 
is regular all the way to the center and reveals a small core that resembles a black hole. 
Further work is expected to produce an equation of state that can account for the high 
gravitational fields observed at about l 16 cm from the center. The presence of a non 
vanishing mass density appears to prevent the formation of black holes but to permit 
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very massive accumulations of dark matter. (The equations of state used in this example 
are not applicable to standard, visible matter.) 
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